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ADSl-RACr 

The emphasis on increased aircraCt and propulsion control systan integration and piloteti 
slnwlation luis created a need for Ixigher fidelity real tliiw dynamic prqxulsion models. A 
real tlnve propulsion systan nodding technique wliich satisfies this need and wliidi provides 
die capdiilities needed to evaluate propulsion systan performance and aircraft systan 
Interaction on manned flight sinulators lias been developed at NASA-Lcwis and daionstrated 
using fli^t simulator facilities at NASA-Anjos, 

A piecewise linear state variable tedinique is used, lliia tedinique provides difi systent 
accuracy, stability and transient response required for integrated aircraft and proiwlsion 
control system studies, Tlie real time dynamic nxxlel laicludea die detail and flexibility 
required for die evaluation of critical control jiarameters and propulsion cai^xnent limits 
over a limited flight enveloiie. Tlie nxxlel contains apprcadmatoly 7. OK Ijytoa of in-line 
caqxitational code and 14. 7K bytes of block data. It tuts an 8.9 tns cycle tino on a Xerox 
Sigma 9 cati{>uter. 

A Pegasus-Harrier propulsion systan was used as a baseline for developing die madiaaatical 
modeling and slmilation tedinique. A hyclroniedianical and water injection control systan was 
also slnulated. Ilie nxxlel lias been programied for interfacing v^ridl a Harrier aircraft 
simulation at NASA-Anes, 

Descriptions of die real time mediodology and nxidel capabilities are presented. 


iNTRODUCriON 

Simulation, widi its tnlierent flexibility, will play a key role in die developnxait of 
Integrated aircraft-propulsion control systans. 'lliese sinulations will provide a 
caiprdiensive source of qualitative and quantitative information regarding the 
diaracteristics of aircraft and propulsion ^stoim in a dynamic state, lliey will also serve 
as tools for die analysis and ayn diesis of control logic and adt test vehicles for control 
software and liardwaro development. 

Since die advent of piloted alaulators and die growing aaphasis for systans integration, 
there lias been an increasing need for higher fidelity real-time propulsion systan nxxlels. 
Propulsion and integrated pontrol ays ton evaluation of aircraft on flight similators require 
diat propulsion systan sinulations be realistic and include significant dynamics as well as 
In^rtant engine parameters, Such sinulations provide die capability to evaluate propulsion 
systans and dielr interaction widi aircraft controls. One sucli modeling technique to 
accanplish dils aipsMllty is described in (1) and (2) . A efynande digital real-time model of 
an advanced propulsion systan based on a piecewise linear mediodology to meet die san» 
requiranents is described in tliis paper. iTils nxxlel provides die engine-control systan 
accuracy, stability and transient response required for studies whidi might include die 
evaluation of critical control parameters, systan response, systan envirormental effects and 
critical propulsion component aerodynamic, medianical and thraixxiiyntimic limits. Tlie model 
may also be used to analyze propulsion control failure modes and effects . 

A Pegasus 11 propulsion system provided the baseline engine for developing die madioiatical 
tVDcleling iind simulation tedinique. The engine nxxlel is a piecewise linear state variable 
representation which was deriv^ fean a detailed aerothermodynamlc simulation of a typical 
Pegasus 11 engine, l^^namics included in. die sinulatlon are engine fan and compressor rotor 


dynamics, engine bumer hcac transfer dynandca aiul engine cwitrol dynwilca, 'Die model 
calculates tranaiejit performance by mnerical integration of t±i>e-dependent differential 
state equations jind contains tlie dynamics necessary to sbrulato aircraft forces rcsultljig 
frcni engine thrust levels, llrLs hlglier fidelity proixilsion control iiimlel provides steady 
state and transient cliaracterlstica for various cngh^c pressures, temperatures, flows, stall 
margins and tltrust. 

Application of the notlol to a sinulatcd flight program and evaluation results are presented 
in (3). 


R12AL ITME MlimiODOUXJY 

Itie real tlaie modiodology is based on a piecewise linear state viiriable techniquo (4) . 
Witliin tliils process, tlae engine nxxlel uses state variables and iivitrix^ fomulations to 
represent tlie engine process at sjxiciflc operating points, 'llie effort presentcxl here was 
directed toward applying tills meUwdology to ncdel gross enghw transients accurately and 
efficiently. Ihe metliodology includes Hie following stepsj 

(1) apply a modal iwalysia and sensitivity study of ulio detailed 
base model to select nxxlel states and input controls 

(2) optimize state variable nxxlel selection to accurately define 
steady state and transient diaracteristics over a specified range 

(3) generate accurate nxxlel pjurtial derivative matrices 
using offset derivative tedmiciuos 

(4) connect and sduxlule die state variable nxxlel 
partial derivatives to represent gross transients 

(5) oiploy a sinple but accurate integration sdxsne for fast computation 

(6) optimize progranming for real time operation 
State Variable Representation 

The state variable representation is sliown in figure 1. It is diaracterized by die follcwlng 
two equations; 

X “ A X + B U (1) 

Y ■» C X + D U (2) 

Equation 1 is a linear constant coefficient natrix differential equatim diat represents 
conputation of engine dynamics. Equation 2 is a linear algebraic equation diat represents 
computation of die observed engine parameters. X is die vector of state variables, X is' the 
time derivative of die state variables, U is die control input vector and Y is die vector of 
observed or output parameters, A is die plant natrix. Its elanenta are die partial 
derivatives of die tina derivatives of each state variable to eadi state variable. Elements 
of the output natrix C define die effect of eadi state variable on eadi outixit variable. Tlie 
control natrix B mid die direct couple matrix D define the effect of eiich control variable on 
eadi state variable time derivative and eadi output paranater. 

Modal Analysis 

A modal analysis is used to determine wliidi states in die nonlinear aerodieimodynamic model 
would adequately roprc.sent die oystwn within the required control bandwiddi. The nonlinear 
aerodiernxidynamic model is dien linearized to obtain die system A, B, C aixl D naixices, An 
eigenvector-eigenvalue analysis of the A matrix is performed. The eigenvectors are examined 
to associate eigenvalues widi states. Higli frequency states outside the control bandwiddi os 
sliown (XI figure 2 are eliminated, A mode concrolldjility matrix vdiidi defines the effect of 
oontrol inputs on states was also generated. States wliidi are uncontnrollable by die inputs 
are eliminated. Hie result is a set of state variable vectors for die real time model. 

hkxlel Resolution 

Modeling gross transir^nt excursions efficiently and accurately in die state variable form 
depends on the nurber of nxjdels selected. Initially, a piecewise linear fit of die steady 
state operating line is performed to define a minlnun resolution. These ttxxlels are dien 
augmented with additional models to accurately define transient response dirougli die full 
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power range and aty excranely nwlijiear areas, 
Matrix ttotlal Derivative Generation 


State variable tedaiiqucs can bo employed to obtain a linear approximation of a nonlinear 
syston by considering operation in die vicinity of a particular operating jxilnt, Matrices A, 
B, C and D ore clwacterized bys 

“IJ "-Sig 

(3) 

axj 

'^ij " 5Uj‘ 

(A) 

Oy. 

"Oxj 

(5) 

auj 

(6) 




To generate finite differoico approxlimtions of tlio partial derivatives, die steady state 
level of each state variable and system input ia indepeixiently stepped in bodi a positive and 
negative direction vrtiile bolding die odier state variables and system inputs fixed. 
Corresponding i»aluos of die state variable derivatives and system outixits are recorded for 
eadi step change. Each matrix partial derivative is calculated by taldng an average of bodi 
step values t 
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This process was autanated on die detailed nonlinear base nxxlol. First, eadi X is perturbed 
one at a time while holding all other states and control inputs constant, 'iliis allciws 
calculation of the A and G matrix partial derivatives, 12adi U ia dien perturbed one at a 
time while holding all odier controls and states constant. Ihis allcws calculation of die B 
and D matrix partial derivatives . 

Several differait levels of perturbations on die states mid Inixits iure used. Perturbation 
step size was minimized to prevent driving model parameters out of range but made large 
enough to excite each state variable tuid output ixirameter. Tlie best overall agreement 
botwemi die linear base nodel and die nonlinear models occurs wlien die states are perturbed 
about 0.5 percent and die inputs about 3.0 percent. 'Ihe partial derivatives generated by die 
offset derivative teclinique are reasonably accurate for steady state at eacli condition for 
wliidi diey were generated, For large transients, ha,;ever, diis is not necessarily true and a 
forced steady state natcli may be required to enswe steady state accuracy fran one state i 
variable model to Uio next along the engine oiierating line. To avoid this, all perturbations 
ore measured fran a steady stote operating lino n^ol wliidi insures .steady state accuracy 
over d.o wliole range. 

State Scheduling Parameter 

llie state variable models must be cormccted efficiently to provide continuous operation over 
die entire range of die ncdel. 'Hie interpolation is contarolled by sdieduling die matrix 
elensnCs widi mi independent variable called the state scheduling parameter (SSP) ! 

SSP « f (X) (8) 


This parameter is cterived from the state variable vector and indicates die relative energy 
level of die engine system. Application of diis procedure results in reducing several linear 
models to one nonlinear model as shewn in figure 3. 


State Veetcar Integration 

A simiile self starting and reasonably accurate integration sdieme is required to minimize 
excessive function evaluation and computation time, Euler integration was found to meet 
diesG requirements and was subsequently a^iopted for diis model development. 
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PRDGRAM DESCRIPTION 


Figure 4 contains an overview block diagram of the inportant state scheduled parameter state 
variable model logic. 

The code contains first pass checks and initialization steps. The initial steacfy state point 
is calculated frcm the state and output operating lines, llie Initial time point of ary 
txansient run is assumed to be in a steady state condition, which requires that: 

AAX + B4U - 0 (9) 

Since the state derivatives will not be zero on the first pass if the initial point is not on 
the operating line, nultiple passes are made through the state dynamic calculation loop until 
the states settle, The solution will be a set of steady state values wliere tiie AX vector is 
given by: 


AX - A"^ BAU 


( 10 ) 


Transient operation occurs as follows. The last time step values for the states are used to 
calculate the SSP. Bracketing model points are then identified relative to the BSP level. 
The SSP is also used to schedule the A, B, C and D matrix elements. These matrix elements 
are stored in a linear equation form (point-slope) which allows for rapid caiputation between 
the discrete model points, AX' a and AU's from the model points above and belcw tlie SSP are 
ccnputed. A relative distance wei^tlng schone is used to carbine deltas from the model 
point above with deltas to the model point below. This basepoint smoothing is set by the 
model point distance from the currait SSP value. 

State derivative carputations are performed by the matrix multiplication of the A and B 
elements with the AX's and AU's conputed earlier. The derivative vector is then Euler 
Integrated and sunrned into tlie storage vector holding the cumulative state level. The 
cumulative AX is . passed out of the matrix integrator and vector sunrned with the basepoint 
values. This final result is then the actual state level to be used for both output and 
calculating the new SSP. 

The basic values for the output parameters are set through the operating line curv'es as 
fijnctions of the SSP. Base output levels change as the SSP dynamically varies investing 
engine energy level changes. In addition to the basic transient response, deviations from 
the output operating line must be calculated, This is done through the output equation which 
xjses the same AX and AU vectors used in the state derivative equation, Tlie aY vector 
computed then represents the change in output values from tlie current set of operating line 
output values . Sunming this value and AY elanents gives the required total output response 
vector. 


APPLICATION 

The methodology presented here was applied to the Rolls-Royce Pegasus 11 propulsion system 
(3) . This system is used in the Harrier VSTX3L aircraft. The primary purpose was to generate 
a real time digital model for use at the NASA-Ames flight sinulator facility to evalu-^te 
integrated control concepts for VSTOL aircraft. 

The source for the real tin:e model development is a con^irehensive aerothermodynaraic 
simulation of the Pegasus 11, Data from Rolls-Royce was used to represent ndnisum engine 
steady state performance on all parameters. A band of 3 percent for specific fuel 
consunption and 1 percent for tlirust and jet pipe tenperature was established as an 
acceptable matching tolerance. Four different match points were diosen covering the full 
power range. As shown in figure 5, the results obtained at these points on the 

aerothermodynamic sinulation ware within 1 percent of the minimum Pegasus 11 steady state 
performance. 

Transient sinu? stion accuracy was obtained by adjusting rotor inertias to specific Pegasus 11 
values, Tlie results are shewn in figure 6. The apparent differences in the initial steady 
state are attributed to differences in the steady state and transient Pegasus 11 data. The 
transient accuracy was coiputed by conparing the maximum rates of change. This was within 5 
percent for all engine parameters. 

Tlie state variable input ar«i output vectors used in the real time model are shown in 
figure 7 . The state vector was derived from the modal analysis , The iiput or control vector 
and output vectors were determined primarily by fuel control and flight envelope 
requiranents. Fourteen linear models are used to define the Pegasus 11 operating 
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charocterlacicB frcm grtxmd idle ac 7 perecsat: power to full uvixiiiw at 109 percent power, 

TI\obo point nodela are linked togetlter by sdiedaling tlie ttatrlx ulcmcaita at each model as a 
functlcan of tlw SSP, In Uiia tnodel the SSI’ waa detinod oa tJvs average of tlwi ateady atate 
fuel flows tliat would occur at oadi jmrticulor current 8 tote viiluaj 



l»l 


Huifl, ateady state fuel flow beconvea an intermediary state sdvodullng parameter derived from 
tl^Q state vector. 


Ibo resulting fliglit sinulator state vartoble engine uodel is a high fidelity proix»l8io«\ 
model whidi provides Steady state mvd traiviient dmracteriotics for desired engine pressures * 
tai^xnraturos, airflows, surge nvrrglns, tlunsts and rotor speeds, 'Ibe canHiuor program 
slnulatlng die engine ailculates botli ateady state and dyiv#uic engine diaructoristics tint 
^^ra ropr<jaenuative of tlio Pc8as>4a 11 engine, 


RESULIS 

Ttie operational real tine digital model contains appreodmatoly 7. OK bytes of In-line 
canputational code and 1A.7K bytes of blodt data, I’rogrannling tedmiciuea doveloiied in (1) 
and (2) were used to optimize Uie cale for real time ojxsration. ibo block data contains all 
of die ABCD matriit daui for die V\ models and aupjxirting functimis . Wttli minor dmnges in 
die code, die nuiber of models used could easily bo dianf’Od. 

For a stoulation tbiie sftep size of 50ms the propulsion system model executed iii 2,3ms on a 
Unlvac 11X0 for a real-to-exocutlon time ratio of 21,7. On die Xerox Sigma 9 canxitor used 
at NASA-Ames, die execurion time was 8.9ms for a ratio of 5,6. Tliese caiputauion times 
Included a detailed hydromeclumical control and mi aircraft force and lialnnce section. 11V3 
engbic portion consuneti about 75 percent of die total execution time, 'llio matrix operations 
iriwlved in die state variable ecpiatlotiS cenSvimeu AO percsit of diat tlim* . 

A lumiber of sots of Crmislent test cases were run to da«onsta.‘nte the real time digital 
Pegasus 11 mcxlol accuracy and capabilities, Only die engine section was osed here, Ibe test 
cases show response to full range 3 second fuel flov'j ron\i transients, 'Ibe minbiMii to maxinugii 
fuel flow rmip provides a moderately rapid transient disturbance diat allowa the full ;et of 
engine partial derivatives to bo exercitsed, A test case at eadi come? of die fll^t 
taivolope as sliown in figure 8 is considered. An additional tijst mise at sea level static 
witli response to the simie fuel floi^ input but including a water injection step at 92 percent 
fan speed was also run. lata fran eadi case was au^mred widi die aerodiermtxlyivimie 
ainulation fran \dvidi it was derived. All parameters {ire in percent of full rmige. 

'J'ransients are si ^n In figures 9 Ou<l 10 for only die full rfinge acceleration at 5000 feet 
{iltitude and tlio same case widi wiiter injection at sea level static. 

Tlie full range acceleration altitude run, figure 9, demonstrates die effects of altitude on 
die nodol. li-misitait bdiuvior and liigli power steady state levels lurc very closely miitdied 
especially for the output dims ts. Since tlimst is die prbuary output to die aircraft 
sbiuliition, it is essential to liave dnust represaitative of die real system. Also, uldiough 
not sliown here, differences up to G percent I'^erc csdiiblted in sane Intenvil engine 
poraneters, in particular, jot ptiie taiperature and fan mid caiprossor surge margin, liiose 
differences could bo significiint if dicy were used as a sensed parameter in a control systan. 

Ibo test cose widi water ii\jection, figvaro 10, sliowod mudi die sane acci.a'acy. Tliem was no 
particular difficulty in modeling die steady state or rtynanic effects widi water iivJeccion, 

Results of die odier tost cases, although not shown here, were similar. Hie effect of 
caiblned altitude and speed showed dvlt altitude Vuid die most effect. Rill rmige transient 
imtdving was excellent here witii most differences occurlng again at low power levels . Speed 
effects alone hud a relatively insignificant effect on die state variable model. High mid 
low s teach' state level mid transient profiles imtdiqd extraiioly well. 

Tlie engine model prosmited here was combined widi a detailed representation of a 

hydraneclianical fuel control mid integrated into a siiiulation for die VUurrler aircraft, Tlie 
results of that evaluation are given in (A) . 
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OONCLUOING RaiARKS 


A stuice scheckiled state variable nKtJwdol<^ has been used fo generate a real time d%ltaX 
prcipulsiw systan sinulatlon for piloted smilators. Itie neUwdology proved effective in 
gavjrating a model of excellent accuracy owr a limited operational range as comiwred to tlio 
oerorhermodynanic model from Wlilch it \*is derived, 'llw inotbodolog/ providea a very floxiblQ 
means to acconpllsh a real time model in a reasonably sliort time, 

llw model executed 5.6 times faster tlum real tln« using an integration time step of 50ms. 
Computiition tin»a analysis lixiicated tluit approxiawtely 30 percent of tlie cycle time \*rs 
consused in matrix operations. It appeaxrs tlmt tlvis could bo a liigl) p^/iyoff arcai for fnturo 
development. Altlioogh tlw total omoutit of code gaioruted in tlxla model is high, 21. 7a bytoa, 
tlvis is generally not a liiiiting factor in real time sbiulntions, 
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Figure 2. - Modeling frequency range of interest. 
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Figure 4. - State scheduled state variable engine modeL 




















CONTROL VEC'IR 

'•FUELaOW 
•WATER INJECTION FLOW 
. • REACTION CONTROL BLEED aOW 
•ENGINE FACE PRESSURE 

• ENGINE FACE DELTA TEMPERATURE 

• AMBIENT PRESSURE 


STATE VECTOR 


AX« 


• LOW ROTOR SPEED 

• HIGH ROTOR SPEED 
•BURNER METAL TEMPERATURE 


OUTPUT VECTOR 


AY ^ 


•CORE NET THRUST 
•FAN NET THRUST 

• FAN PHYSICAL 

• AIRFLOW 

• HIGH COMPRESSOR PHYSICAL AIRFLOW 

• FAN AVERAGE PRESSURE RATIO 

• HIGH COMPRESSOR PRESSURE RATIO 

• FAN NOZZLE DUCT TEMPERATURE 

• FAN NOZZLE DUCT PRESSURE 

• HIGH COMPRESSOR INLET TEMPERATURE 

• HIGH COMPRESSOR INLET PRESSURE 

• BURNER INLET TEMPERATURE 

• BURNER INLET PRESSURE 

• HIGH TURBINE INLET TEMPERATURE 

• HIGH TURBINE INLET PRESSURE 

• LOW TURBINE INLET TEMPERATURE 

• LOW TURBINE INLti PRESSURE 

• CORE NOZZLE DUCT TEMPERATURE 
•CORE NOZZLE DUCT PRESSURE 

• THRUST SPECIFIC FUEL CONSUMPTION 

• FAN SURGE MARGIN" 

• CORE SURGE MARGIN" 


«-AUXILLARY MODEL OUTPUT COMPUTATIONS 


Figure 7. - Real time model state variable vectors. 



Figure 8. - Real time model flight envelope. 
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